Mon. Not. R. Astron. Soc. 000. [THTTl (2003) Printed 2 February 2008 (MN IATeX style file v2.2) 



Cross Terms and Weak Frequency Dependent Signals in 
the CMB Sky 

C. Hernandez-Monteagudo 1 * and R.A.Sunyaev 1 ' 2 

1 Max-Planck Institut fur Astrophysik, Karl-Schwarzs child- Str. 1, D-85740 Garching, Germany 

2 Space Research Institute (IKI), Profsoyuznaya 84/32, Moscow 117810, Russia 



ABSTRACT 

In this paper, we study the amplification of weak frequency dependent signals in 
the CMB sky due to their cross correlation to intrinsic anisotropies. In particular, we 
center our attention on mechanisms generating some weak signal, of peculiar spectral 
behaviour, such as resonant scattering in ionic, atomic or molecular lines, thermal SZ 
effect or cxtragalactic foreground emissions, whose typical amplitude (denoted by e) 
is sufficiently smaller than the intrinsic CMB fluctuations. We find that all these ef- 
fects involve either the autocorrelation of anisotropies generated during recombination 
(z rec ) or the cross-correlation of those anisotropies with fluctuations arising at redshift 
Zi. The former case accounts for the slight blurring of original anisotropies generated 
in the last scattering surface, and shows up in the small angular scale (high multi- 
pole) range. The latter term describes, instead, the generation of new anisotropies, 
and is non-zero only if fluctuations generated at redshifts z rec , Zi are correlated. The 
degree of this correlation can be computed under the assumption that density fluctua- 
tions were generated as standard inflationary models dictate and that they evolved in 
time according to linear theory. In case that the weak signal is frequency dependent, 
(i.e., the spectral dependence of the secondary anisotropies is distinct from that of 
the CMB), we show that, by substracting power spectra at different frequencies, it 
is possible to avoid the limit associated to Cosmic Variance and unveil weaker terms 
linear in e. We find that the correlation term shows a different spectral dependence 
than the squared (oc e 2 ) term considered usually, making its extraction particularly 
straightforward for the thermal SZ effect. Furthermore, we find that in most cases 
the correlation terms are particularly relevant at low multipoles due to the ISW effect 
and must be taken into account when characterising the power spectrum associated 
to weak signals in the large angular scales. 

Key words: cosmic microwave background - large scale structure of Universe - 
galaxies: clusters: general - methods: statistical 



1 INTRODUCTION 

Standard theories state that the field of density per- 
turbations arising after the inflationary epoch, (<5(x) = 
(p(x) — p)/p, with p the average density), should be gaus- 
sian, homogeneous and isotropic, jGhith|^)8l|; IStarobinskiil 
ll98ll:lMukhanov fc Chibisovll982tlLindell983l) . The Fourier 
modes of this field are predicted to have indepen- 

dent real and imaginary components, which should be 
gaussian distributed from a scale- invariant power spec- 
trum, (Harrison-Zel'dovich, (HS), (IZeldovicbl Il972f) ^. i.e., 
{S k Sq) = (27r) 3 P(fc)fo(k + q), with P(k) oc k. This power 
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spectrum determines the properties of the spatial correlation 
of the perturbation field, since it is the mere Fourier trans- 
form of the correlation function. These perturbations are 
small compared to the homogeneous background, (\S\ <C 1), 
but grow up due to gravitational instabilities. This growth 
is independent for each mode, i.e., mode coupling can be ne- 
glected, as long as the perturbations remain small and linear 
theory can be applied. 

From the observational point of view, the first test 
ground for this perturbation field is the study of the tem- 
perature anisotropies of the Cosmic Microwave Background 
(CMB). Most of these temperature fluctuations were gen- 
erated by energy density inhomogeneities in the universe 
during the epoch at which most electrons and protons re- 
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combined to form hydrogen and radiation decoupled from 
matter, (last scattering surface, LSS). At this stage, the 
density inhomogeneities were still under linear regime, pro- 
vided that the amplitude of typical measured CMB temper- 
ature fructuatjons are o n e part in one hundred thousand, 
(e.g. ISmoot et all (Il99lf) : iBennett etaH (l2003al) '). In their 
transit from the LSS towards us, the CMB photons wit- 
nessed the matter collapse and formation of non linear struc- 
tures such as galaxies, clusters of galaxies, filaments and 
superclusters of galaxies that today conform the visible uni- 
verse. The crossing through these scenarios imprinted on the 
CMB photons new temperature anisotropies, which are usu- 
ally labelled as secondary. The amplitude of these secondary 
anisotropies is, in many cases, few orders of magnitude be- 
low the level of the primary ones generated in the LSS. How- 
ever, the different spectral behaviour of some of them might 
help in the distinction from the primary. In this context, the 
presence of foregrounds, galactic or extragalactic, with their 
own spectral dependence, will make the picture furtherly 
more complicated. 

Consequently, a major issue in current CMB science 
is the accurate component separation in future microwave 
maps. From the observation point of view, a set of space 
and groundbased experiments with unprecendented sensi- 
tivity and angular resolution, counting with several broad 
band detectors spread in appropiate frequency ranges, are be- 
ing p roposed or alread y under development, (e.g., Planck 1 , 
ACT jKosowskvl 120031) . South Pole Telescope 2 , QUIET 3 or 
CMBPOL 4 ). In the theoretical side, new analysis techniques 
of temperature maps, based both in real and Fourier space, 
and dealing with second (correlation function, power spec- 
trum) and higher order momenta of quantities derived from 
the temperature field are being developed and tested on sim- 
ulated data. Nevertheless, the two main limiting factors in 
this task will be i) the instrumental noise and instrument 
systematics and ii) the cosmic variance, associated to the 
fact that our characterization of the universe is statistical, 
but based on a single realization of it. 

In this work, we study the spatial correlation of density 
fluctuations in the universe, and how this reflects in the 
CMB angular power spectrum. These aspects must be 
taken into account if an accurate characterization of the 
CMB power spectrum is to be achieved, particularly at the 
large angular (low multipole) scales. This study also allows 
us to propose a method that uses observations in different 
frequencies and combines power spectra in such a way 
that avoids the limitation imposed by cosmic variance, and 
unveils weak signals whose amplitude a w is in the range 
o~n > cr-ui > o~m (crjv /o~t), where crjv is the experimental 
noise amplitude and a t is the typical amplitude of a 
dominant signal t which is assumed to be totally correlated 
to the weak signal w. This approach was al r eady utilized 
in iBasu. Hernandez-Monteaeudo fc Sunvaevl i2004l) when 
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charactering the effect of metal atoms and ions on the CMB 
during the secondary ionization. 

In Section 2 we outline our method, which we apply 
in Section 3 on particular physical mechanisms generat- 
ing secondary fluctuations in two different cosmological 
scenarios. In Section 4 we comment our results and conclude. 



2 COMPARING SECOND ORDER MOMENTA 
2.1 The Flat Case 

Our starting point will be the superposition of two signals, 
gi(y) t\ and g^iy) ti, whose amplitudes will show, a priori, 
a different frequency (v) dependence. In real space they give 
rise to 



T(x,i/) = 0i(i/)ii(x) +32(f)e*2(x) +N„(x), 
and analogously, in Fourier space, to 

T kM ^iH^.k + SaHrfak + ^k' 



(1) 



(2) 



where k is the Fourier mode under consideration. If we as- 
sume that <7i(^i)ii and giiy^tz are of similar amplitude, 
then the parameter e gives the relative amplitude of both 
signals , and for the cases considered below, we shall take 
e«1. N v is the noise component present in the map. Now 
let us assume that the experiment is able to observe at two 
different frequencies V\,V2- Defining / = gi{fi) / '31(^2) , we 
find that: 

<5(T 2 (x,y))= T(x,^)T(y,^)-/ 2 T(x,^)T(y,^) = 

giiyi) gi(vi) 



e 31(^1)32(^2) 



+ e 2 gl(yi) 



g\(vi) 32(^1) 

giiyx) 32(^2) 
01(1*2) ffa(fi) 



[*i(x)t 2 (y) -Ki(y)* 2 (x)] 
fe»(x)ta(y) 



+ N V1 (x)N vl (y) - f N V2 (x)iV, 2 (y) 
+ 0[N V1 ,N V2 ], 

or, in Fourier space, 

8 (T£ q ) = T k MTqM - / 2 T k (/, 2 )T q (/, 2 ) = 



(3) 



e £1(^1)32(^1 
+ e 2 32(^1) 



3i(^i) 32(^2 
32(^2) 32(^1) 

31(^1) 32(^2 
32(^2) 32(^1) 



[t^kfe-q + * a ,k*i.q] 

*2,k*2,q 



+ ^ 1 ,k^ 1 .q-/^ 2 ,k^ 2 .q 

+ 0[N V1 ,N V2 \. 



(4) 



O [N VI , N U2 ] in both equations refers to cross terms of the 
noise field with all the other components at a given fre- 
quency. These two equations should be compared to the 

squared difference map, ((ST (x, y)) 2 , f<5?kq^ )> gi ven by: 
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(<5T(x,y)) 2 ^ 

(T(x, Vl ) - /T(y, v 2 )) (T(x, Vl ) - fT(y, u 2 )) 
e 2 ( [g 2 { Vl ) t 2 (x)] 2 +[/ g 2 (v 2 ) t 2 (y)] 2 - 

2/52^l)ff2(^2)t2(x)t 2 (y) ) + 

7V 2 1 (x)+iV 2 1 (y)+0[7V I , 1 ,iV, 2 ], 
in real space, and 

5T, 



(6) 



(5) 



(T k (^) - /Tq(^)) (T k (^) - /TqM) 

e 2 (^2,xff2(^i)] 2 + [/t 2 ,yS2(^2)] 2 - 

2 / 32(^1)32(^2)^2, xi 2 ,y ) + 

in Fourier space. 

It is clear that for e < 1, 8 (T 2 (x, y)) or 8 (t^ q J are 

much more sensitive to the weak signal et 2 than (<5T (x, y)) 2 

or (^ST^ . The obvious difference is the term linear in e 

present in eas. l|3l4p . However, in the context of Cosmology 
and CMB, one counts with only one single realization of the 
Universe, and the quantities defined above as 5 (T 2 (x, y)) 

or 5 (^Ty, must be averaged either in real or Fourier 

space, in order to acquire some statistical meaning, (i.e., if 
averaged under certain conditions, they yield estimates of 
the correlation function and the power spectrum, respec- 
tively). After this average, the term linear in e becomes 
proportional to {t\ t 2 ), and will not average out if and only 
if both signals ti, t 2 are correlated, at least to some extent. 
Therefore, in order for this cross term to be of any utility, 
both the dominant and the weak signals must be correlated. 
We shall show below that this is indeed the case for signals 
coupled to linear fluctuations of the density field generated 
after inflation. Another point to remark is that, because of 
substracting quantities computed from the same maps, one 
exactly cancels the dominant signal, leaving no room for 
the uncertainty due to the cosmic variance associated to it. 
This allows the weak signal be under the limit imposed by 
the cosmic variance of the dominant one. 

As mentioned in the Introduction, in linear theory 
all Fourier modes 5-^ of the density fluctuations evolve 
independently according to a growth factor D(n) (77 is 
conformal time) which is dependent on the cosmological 
parameters of our universe. These modes are all indepen- 
dent, and for reasons associated to the homogeneity and 
isotropy, must depend exclusively on the modulus of the 
k vectors, k. This allows writing the power spectrum as 
(^(Jq) = (27r) 3 P(fc)5_D(k + q). In an analogous way, the 
averages of the product of all pair of quantities depending 
linearly on 8^ will be proportional to the power spectrum. 
This applies practically to all perturbations of physical 
quantities, such as peculiar velocities or gravitational poten- 
tials, that are responsible for the generation of temperature 
anisotropies in the CMB. 



The average in our maps will be performed in the real 
space in such a way that the distance between x and y is 
kept constant. In Fourier space, we shall take 5 q equal to 
k, fix the modulus (k) and average over the mode phases. 
The former will yield the correlation function, the latter the 
power spectrum. This average also removes all cross terms 
in noise. Furthermore, if we assume that the statistical prop- 
erties of noise have been characterized, then it is possible to 
substract the expectations for the terms quadratic in noise 
in eas. 131411 . and the residuals of this substraction can be 
treated as random variables. These random residuals should 
be regarded as the effective noise in our correlation func- 
tion or power spectrum estimates, and will be denoted by 
AA^.x-y and AN v , k : 

AA^.x-y = £((N v (x)N v (y))) - (N v (x)N„{y)) E xp (7) 

AN„, k = £ ((iV^.k)) - (N^N^bxp, (8) 

where £ and the label EXP denote estimated on the map and 
expected values, respectively. We are assuming that noise in 
different frequencies is uncorrelated. For the case of gaussian 
white noise, it is easy to prove that (AA^x-y) = {AN v ,k) = 
0, (assuming a correct characterization of noise), and that 



<A< X -y) = -(N u {x)N v (y)) 



<A7V 2 fc > = -(JV„(k)JV„(-k))! XP . 



(9) 
(10) 



n is the number of points, either in real of Fourier space, 
used when estimating the averages 6 . Having this in mind, 
we can perform the averages and rewrite eas. 13141 like 

£((S (T 2 (x,y))» = 

gi(vi) gn{vi) 



g\(v x )g 2 (v 2 ) 



± A/ 



and 
£ 



H?k,-k)» 
£31(^1)32(^1) 



± A, 



31(^2) 32(^1) 
+ O [e 2 ] 



3i(^i) 32(^2) 
31(^2) 32(^1) 

, + O [e 2 ] . 



2 (*i(x)ta(y)) 



(11) 



(*x,k*2,-k + *2,k*i,-k) 



(12) 



Ajv is the residual noise contribution, A^ = {AN 2 ^ + 
f 2 (AN 2 2 ) in both real and Fourier space, (eqs.llS HTH '). From 
this equation, one can see that the approach proposed here 
will be sensitive to e t 2 if: 



a N 

lu e > — X 



yJl (! + / 2 ) 



21 



91(^2)92(1^ 



21 V 
J) 



(13) 



with aff = (N 2 (x)} taken equal for the two frequencies and 
of = gf(ui)(t 2 ) ~ 3 2 (^2)(t 2 ). The factor u accounts for the 
cross correlation between t\ and i 2 , i.e., u> = (ti ■ t 2 )/a 2 , 
(note that in the absence of correlation, u) — since we 



5 Note that, for real signals, <5_jj = <5j^. 

t noise, if x 
{(^(x)) 2 )((^(y)) 2 >/n 



6 For white gaussian noise, if x ^ y, then (A7V 2 x _y) 
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are taking (ti) — by construction). Let us remark that the 
limit on ue is roughly an order in ajv/crt beyond the limit im- 
posed on e by eqs. H5I6|I . Note that in the case of similar fre- 
quency dependence for the two signals, (gi(u) ~ Q2{v)), this 
method cannot work. For similar reasons, if g\(v) 7^ (72 {f), 
then it should be possible, a priori, to perform as many 
consistency checks in different frequencies as the instrument 
permits, since the correlation term should vary its amplitude 
as dictated by the frequency dependent term 



SCi = 2eY{y 1 ,v 2 ){ai, m {af^ k Y) + o[t 2 



(17) 



r(i/i,i/ 2 ) = 31(^1)32(^1) 



1 



giiyx) 52(^2) 



g\{v 2 ) g»{vi) 



(14) 



From this formalism, it follows that the importance 
of this approach relies i) on the amplitude of the cross- 
correlation between the signals under consideration and 
ii) on their spectral dependence. Let us remark as well 
that this method is sensitive to the relative sign of the two 
signals. In the context of the CMB, this correlation will 
preferrably show up in the low multipole range: at these 
large angular scales the instrumental sensitivity performs 
best, but the removal of galactic foregrounds becomes 
particularly difficult. In the next section, we shall address 
several scenarios where this correlation may be relevant, 
and discuss under which conditions the method proposed 
here becomes useful. 

The approach outlined here is comple- 
mentary, bu t di fferen t, to that used in, e.j 



Bandavetail (11996ft. iKneissl et alj Jl997tl. 

iubin o-Martm. Atrio-Barandela. fc Hernandcz-Montcagudo 
1 200fJ), and more recently, iBoughn fc C rittenden! 

2003|i. |Fosalba 1 Gaztana ga fc C^stan^rl ~~pOO^ and 
[ lernlm^le^^ontea^uck^^^i^^ i2004T) . In all 

those cases, the weak signal (efe) was not considered to be 
correlated to the dominant signal, but it was cross-correlated 
to an external template: this cross-correlation retained only 
the term linear in e, and hence no substraction was required. 

Hereafter, the term proportional to e will be referred to 
as the linear or cross term, whereas the term proportional 
to e 2 will be denoted as the squared term. 



2.2 Correlations Projected on the Sphere 

In this subsection we briefly outline the formalism that de- 
scribes the analysis of temperature fluctuations in the CMB. 
It is customary to work in the spherical Fourier space, in 
which the coefficients a; im 's define a temperature field in 
the celestial sphere through the following decomposition on 
spherical harmonics: 



f(M)-i; 



(15) 



The power spectrum for an arbitrary temperature field 
is obtained after averaging the Fourier coefficients, 



Ci = {ai im ai t m)- 



(16) 



Having this in mind, the analysis of weak signals out- 
lined in the previous section translates into the spherical 
case as 



However, when computing this correlations, it will be 
convenient to express the az, m 's as integrals in the flat 
Fourier space. Indeed, the tem perature field can be de com- 
posed in Fourier modes as fe.g. lHu fc Sugivamal ((,1995)): 



A ( k. 11. //,, ) = / dx^(x, n, 770) e lkx 



5^(-i) i (2t + l)^(/x)A,(k,f, ) 



(18) 



with fi — k • h, and n is the pointing vector on the sky 
given by (9, (f>). 770 denotes the conformal time evaluated 
at the present epoch. The last step shows the expansion 
on a Legendre polynomial basis, and assumes implicitely 
that perturbations are a xially symmetric about k, (e.g., 
iMa fc BertschingeJ il995h ~). From this, it is straightforward 
to show that, for x = 0, the a^ m multipoles can be written 



ai, m = H) ! 4tt / dkY,* m (£) A;(k, 77). 



(19) 



In linear theory, A;(k, rj) — A;(fc, rf)ij)i (k), with 
ipi (k) the initial scalar perturbations and {ipi (k) ipi (q)) — 
(fc) (2it) 3 S d (k + q) the initial scalar perturbation power 
spectrum. It turns out that, after integrating the Boltzmann 
equation, the mode A (k, n, 770) can often be written as a 
line-of-sight (L OS) integral of some sources dependent on k 
and rj, SQl,T}), (ISeliak fc Zaldarriagalll996l : 



A(k,n,7 ?0 ) = / d V e lk ^°' v] S^, V ), 



(20) 



where the sources can be related to the velocity, poten- 
tial and/or density perturbation modes. After using the 
Rayleigh expansion for the exponential in equation 1201 , it is 
eas y to show that the multipol es A;(k, 770) can be expressed 



easy to snow mat tnc multipole 
as ISeliak fc Zaldarriagalll99r3l : 



A ; (k,J7o)oc dr/ji[k(rio - r])]S(k,rj). 



(21) 



This is only correct if the source term has no [t, dependence, 
S(k, 77) = S(k,rj). Otherwise the integral along the LOS is 
projected on spherical Bessel functions of different order, 
(i.e. A; is an integral of ji+i and ji-i if S(k, 77) <x fj,). In 
all cases considered here, the sources will be [i independent, 
and eq. (I21H will be used. 

This expresion of A;(k, 770) also allows us to make some 
predictions regarding the multipole range where the cross- 
correlation term (ai iTn (a]"^ k )*} will be relevant. The formal 
way to see this is through the integral defining SC'i: 



(22) 



5Q oc I dkk 2 P^(k) Si(k,Tji) S 2 {k, Vl ) x 

ji(k[Vo - T)i])3i(k[Vo - m\ 



In this equation, we have assumed that the two signals have 
been generated at conformal times 771 and 772 (with 771 > 772). 
For a fixed I, we have that ji(x) ~ 1 if x ~ I. For x <C 1, 



3l{x) 



and ji(x) ~ cos(a; — lir/2 — tt/4)/x if x ^> I. 



From this it is easy to see that, for a fixed I, the spherical 
Bessel functions will be close to unity if k ~ fci = l/Srji, 
k ~ k-i = l/Sr/2 in each case, (Srji = 770 — rji, i — 1,2). In 
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practice, this means that, given that ki > k2, for the k 
range for which ji(k8r]2) is unity kSrji < I, so that, for the 
very low k's (and hence very low I's), ji(k5r/i) will approach 
to zero if kSr/i <C 1. This reflects the fact that such modes 
do not enter in the angular scales given by I, and it is 
easy to show that this will take place predominantly in 
multipoles below l min = (770 — ip)/(r]o — On the other 
hand, for the k range for which ji(k8r)i) ~ 1, we then 
have that ji(kSrj2) ~ cos(k8ri2 — Itt/2 — tt /4) / (kSrfe) if 
kSrj2 S> I- Hence, the phase difference between both Bessel 
functions will become important if k{r]i — 7/2) ~ 27r, or 
equivalently, for l max ~ 2^(770 - rja)/(»7i - 772)- lmax stands 
for the multipole at which we expect a change in the 
cross-correlation structure between two relatively nearby 
signals. However, we may find scenarios in which both 
signals are so distant that l m ax ~ 1, and for which this 
analysis cannot be applied. Also, we must keep in mind the 
caveat that we are ignoring the k dependence of the sources, 
which condition the actual amplitude of the correlation. 



2.3 Frequency Dependence of the Cross Terms 

We next focus on the frequency dependence of the SCi's. 
This method is based upon the assumption that dominant 
and weak signals have different spectral dependence. This 
translates into a frequency dependence of the SCi's given 
by: 

SCi = (32(^1) -52(1^2)) x 
^2e(a iim «r fe )*> + 6 2 (32^i) + <72(^2))(K^ fc | 2 )^,(23) 

where we have taken ti to be the primordial CMB fluc- 
tuations and hence gi(u) = 1 = const. This equation 
shows the frequency dependence of the SCi's and also 
manifests the different behaviour of the correlation term 
and the squared term with respect to v. That is, if we define 
A(ff) = 32(1*1) - 32(^2) and A(g 2 ) = gl(m) ~ g\(vi), then 
the (linear) cross term is proportional to A(g), whereas the 
squared term is multiplied by A(<? 2 ), e.g., the latter is more 
sensitive to big changes in g{y). This different behaviour 
should motivate the choice of observing frequencies in order 
to distinguish the contribution of both terms. 



2.4 Relative Sign Dependence of Weak and 
Dominant Signals 

Since the cross term couples different signals, it is sensitive 
to the relative sign or phase present between them. That 
is, it is sensitive to whether both signals are correlated 
or anticorrelated. This sign depends upon the physical 
processes relating both signals and their particular spectral 
dependence, and can be different in different I ranges. 

In the case of the thermal Sunya e v-Zer dovich effect 
(hereafter tSZ, ISunvaev fc Zel'dovichl Jl980t) 1. we shall 
find that, for the low frequencies for which the effect 
decreases the CMB brightness, (v < 218 GHz), the tSZ 
will be anticorrelated to the intrinsic CMB temperature 
fluctuations (caused mainly through the late ISW effect), 



whereas for v > 218 GHz both signals become correlated. 

For resonant scattering, at high i's, we shall see that 
blurring of original CMB anisotropies dominates (SCi < 0), 
whereas at low multipoles generation of new anisotropies 
make SCi > 0. 

These scenarios are addressed in detail in the next 
Section, although we stress that this sensitivity to the 
relative phase/sign of the fluctuations is intrinsic to our 
method, and applies to any pair of signals. 

This relative sign dependence leads to the specific 
(angular) /-dependence of the effects under consideration, 
and both aspects show up combined in the final SCi's. 



3 PARTICULAR CASES AND POSSIBLE 
APLICATIONS 

In the context of CMB, the cross term e(ai >m (af^)*} dis- 
cussed above appears due to different physical processes. In 
what follows, we shall analyse the most relevant in two dif- 
ferent cosmological scenarios: the ACDM model suggested 
by WMAP observations, with cosmological parame- 
ters (n m ,n A ,^b,h,n s ) = (0.248,0.798,0.044,0.72,1.), 
and a critical Einstein-de Sitter Universe with 
(n m ,£lA,Slb,h,n s ) = (0.956,0,0.044,0.72,1.), (here- 
after denoted as SCDM ). The inclusion of SCDM model 
responds to the need of understanding the correlations in 
scenarios with no ISW effect. 

The growth of the Large Structure of the Universe is 
such that it is the small overdensities the first ones to be- 
come non linear and form the first haloes, which, with time, 
merge to form more massive structures. In order to see the 
effect of these haloes on the CMB power spectrum one must 
focus on the typical angular distance between sources. If 
sources are distributed uniformly, then one must take into 
account only the so-called poissonian term, but if sources 
are in some way clustered, then a correlation term must be 
also c onsidered, llLacev fe ColdlT993t iKomatsu fc Kitavamal 
1999). These two contributions conform what we have called 
the squared term, proportional to e 2 . 

The approach proposed here provides an additional 
way to study the effect of the halo population on the 
CMB, consisting in looking at the correlation of their 
spatial distribution with the intrinsic CMB temperature 
anisotropies; i.e., the cross (linear in e) term. This coupling 
responds, in most cases (but not all), to the correlation of 
the density fluctuations field with the gravitational poten- 
tial fluctuation field in a ACDM universe, (ISW effect). The 
particular spectral dependence of the cross term compared 
to the squared term makes it feasible to distinguish between 
them, enabling a separate and independent analysis of the 
halo population. 

We must note that the nature of the correlation 
is independent of the particular physical process, but 
hinges exclusively on the spatial distribution of haloes. 
To model the halo population, we have recurred to the 
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Press-Schechter formalism, dPress fc Schechterll974t) . which 
in general provides a good fit to the outcome of numerical 
simulations, although small correc ti ons to it have bee n 
suggested, JSheth fc Tormenl Il999l : I Jenkins etafl 1200 ll) . 
The latter can be easily implemented in our procedure. 
However, this description of the halo population must 
be accompanied by a proper modelling of the physical 
environment in the haloes, which condition the physical 
phenomena under study, (i.e., the fraction of neutral 
hydrogen in 21 cm emission, the cosmological history of the 
star formation rate in dust emission, the number density of 
radio galaxies versus redshift for radio background studies, 
etc). 



3.1 Thermal SZ Effect and intrinsic CMB 
fluctuations 

The tSZ effect arises as a consequence of the Doppler change 
of frequency of CMB photons due to Thompson scatter- 
ing on fast moving thermal electrons. In this scattering, the 
transfer of energy from the electrons to CMB photons trans- 
lates into a distortion of the Black Body spectrum of the 
CMB radiation. Consequently, the tSZ effect introduces fre- 
quency dependent temperature anisotropies in the Cosmic 
Microwave Background, which, in the non-relativistic limit, 
can be written as an integral of electron pressure along the 
line of sight, 



ST 



, , N k B T e (rj) 
dr) a(l]) —a T n e {ri), 



(24) 



with g(x) = xcoth(x/2) — 4 and x = hv/ksTo the adimen- 
sional frequency in terms of the CMB monopole To. For 
this reason, clusters of galaxies, with their gravitational 
wells filled with hot gas acting as sources of electron 
pressure, constitute the main target of tSZ observations. 
However, diffuse ionized gas, placed in the larger scales 
of superclusters and filaments where still some pressure 
support is provided, should also leave an imprint on the 
CMB spectrum by means of the tSZ effect. However, 
this effect is, for / < 2000, remarkably smaller than the 
intrinsic CMB anisotropies, and this allows us to apply the 
formalism outlined above. 

Recently there has been active discussions about the 
origin of some excess power f ound at I ^ 2000 in ground- 
based CMB experiments, iMason fc CBI CollaboratiorJ 
20011 : iGoldstein et~aH 12003ft . Some groups feond et alJ 



2002) have argued that it can be due to tSZ signal coming 
from unresolved galaxy clusters. Since the power spec- 
trum is a quantity which, a priori, does not retain sign 
information, methods based on the sign of the skewness 
of the probability distribution function of the signal have 
been devoloped in order to discern whether such signal 
comes from negative tSZ clusters or positive point sources, 
jRubino-Martin fc Sunvaevll2003l) . 

In what follows, we show how the frequency dependence 
of the SCi's can be of relevance in this problem . We shall use 
an approach similar to that of ICooravl (1200 lT) to model the 
temperature fluctuations introduced by the population of 



galaxy clusters. The fc-mode of the temperature fluctuation 
field is given by the following LOS integral: 



A(k,?7o) = / dn g{x) e 



ifc/4'70-'?] 



dM f( v ,M) 



M 



1/3 



a T n e (r])T e (M,r])V(M,r]) 



b{M,r{) 



x <5 L 



(25) 



n e (n) is the background average electron number density 
at epoch n, T e (M, r?) is the cluster el ectron temperature 
given by, e.g- lEke. Cole, fc Frenkl lll996l) . and b(M,r)) is the 
halo bias factor. (iMo fc White! Il99fj) . TXM.n) is density 
LOS integral for a f3 — 2/3 model (for the case in which 
the line of sight goes through the center of the cluster), 
V(y M) = 2r c (M,?7)tan" 1 (p) jAtrio-Barandela fc Miicked 
1999) , with r c the cluster core ra dius the same as used by 
Rubi no-Martm fc Sunvaevl <l2003l) , and p is the virial radius 
to core radius ratio, which we have taken to be 10. The 
mass integral multiplying represents the pressure bias 
generated at galaxy clusters, and is characterized by the 
mass function f(n,M) (for which we have used the Press- 
Schechter (PS) formalism): 



/(»J,M) = A /J 



da 



DM 



(26) 



where 8 C is the spherical collapse critical overdensity and 
o~(M,T]) the mass fluctuation field. 

In Figure Q we show our results for the ACDM 
universe. The amplitude of the cluster induced tSZ power 
spectrum (square term evaluated at Rayleigh- Jeans (RJ) 
frequencies) equals that of the intrinsic CMB power 
spectrum at I ~ 2000, and then drops steeply due to the 
lack of very high k modes in our integration, (dashed 
line). Nevertheless, we remark that this approach to 
model the cluster induced signal ob serves the effect of the 
clust er-cluster correlation term, llKomatsu fc Kitavamal 
1999), since its dependence versus I is not C; oc const at 
low multipoles, as it would be expected for the poissonian 
term. Provided that, in this model, cluster induced tSZ 
temperature fluctuations are determined by the matter 
density fluctuation field, its correlation properties are also 
governed by the matter power spectrum. Let us also remark 
that there is no flat approximation here, and hence the 
predictions should apply to the very large scales. In the 
small scales, for which the squared term is dominant, our 
model is comparable to the res u lts of N-body simulations, 
JSpringel. White, fc Hernsquisd 1200 it iKomatsu fc Seliakl 
12001 IZhang. Pen, fc Wand 120021) . who, compared to each 
other, provide relatively similar predictions. 

Our approach aims to describe the interplay between 
the linear and the squared terms, together with their com- 
bined effect. However, we do not intend to provide accurate 
predictions for the amplitude of the tSZ-induced power 
spectrum: this is an open issue subject to be explored via 
hydrodynamical simulations and a better understanding of 
the distribution of galaxy clusters with respect to redshift. 
Progresses at this respect should leave our qualitative 
descriptions of the frequency and / dependence of the tSZ 
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Figure 1. Sign and frequency dependence of tSZ fluctuations. 
(In both panels:) Both dashed and dotted lines give power spec- 
tra in the RJ range, for which the tSZ shows no frequency de- 
pendence. The power spectrum of the squared term associated 
to the tSZ effect due to the cluster population is given by the 
dashed line, (ACDM model). For the sake of comparison, the 
CMB power spectrum given by the cosmological parameters pro- 
vided by WMAP's team is displayed by the upper solid line. 
The dotted line gives the absolute value amplitude of the cross- 
correlation term between the intrinsic CMB and the tSZ signal. In 
the upper panel, the (bottom) solid lines give the actual predicted 
|<5C?;|'s obtained after taking a 218 GHz channel as reference, for 
two close observing frequencies (100 GHz and 143 GHz). We have 
assumed that the tSZ signal cancels exactly in the reference chan- 
nel, and that the SCi's are entirely due to tSZ effect. Note the 
change of sign of the total power (linear term plus squared term) 
at hero, below which the linear term dominates, due to the an- 
ticorrelation of tSZ signal and CMB at low frequencies. The dot- 
dashed and the three dot-dashed lines give the nominal amplitude 
of the noise residuals for the HFI 143 GHz and 100 GHz channels, 
respectively. They are well below the signals we are studying. In 
the bottom panel, we take 230 GHz, 270 GHz and 545 GHz as 
observing frequencies and 218 GHz as reference channel, finding 
no change of sign for the <5C;'s, (bottom solid lines). 



8=4000 



1=40 



i=8 
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100 200 300 400 500 
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Figure 2. Frequency dependence of the SCi's for I = 8,40 and 
4000. For I = 8, SCi behaves versus frequency as gtszix) = 
iccoth(:r/2) — 4, whereas for I = 4000 the v scaling is propor- 
tional to gl sz , and never crosses zero. The I = 40 is a linear 
combination of these two extreme cases, weighted by the relative 
amplitudes of the linear and squared terms. 



power spectrum untouched. 

We can see in figure Q that the absolute value of 
the cross term evaluated at RJ frequencies (dotted line) 
shows an amplitude a factor 5 to 20 higher in the large 
scales (I < 20) than the dashed line (squared term). Once 
the frequency dependence of the cross (linear) and squared 
terms is taken into account, we find different patterns 
for the SCi's according to the observing frequencies. For 
v < 218 GHz, we see in figure Q that the SCi's become 
negative in the low-Z range for which the linear term 
dominates, and the particular multipole at which SCi's 
cross zero (hereafter referred to as l zer o) depends also on 
the observing frequency. The value of such multipole for 
different frequencies in the ACDM model is shown in figure 

it remains roughly constant in the RJ regime, but 
approaches higher values as the frequency tends to 218 GHz. 
This is due to the fact that the squared term tends to zero 
much faster than the linear one when frequencies approach 
218GHz. For v > 218 GHz both linear and squared term 
are positive and hence SCi does not change sign. Note 
that these predictions for the SC'i 's versus I and frequency 
are specific only for the tSZ effect, and should permit to 
distinguish it from the contribution of other sources. 

The different dependence versus v for different l's is 
displayed in Fig.(|5J: the two extreme cases are given for 
SCi at I — 8 and 4000, whereas the intermediate case 
corresponds to I = 40. The behaviour of the SCi's versus 
frequency is a consequence of i) the independence of the 
photon spectrum upon redshift and ii) the fact that the 
tSZ surface brightness changes sign at v = 218 GHz. 

After defining the correl ation coefficient as 
TZi = (af : ^ B aj s ^)/ v /Cf M B Cf z , we plot it for both 



8 C. Hernandez- Monteagudo & R.A.Sunyaev 



. neqative sources 



positive sources 
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Figure 3. Multipolc at which SCi induced by the tSZ ef- 
fect change sigh, for observing frequencies below 218 GHz, 
for which clusters can be regarded as negative point sources, 
jKorolev. Svunvaev. fc Yakubsevl Il986h . Note that it remains 
constant for RJ frequencies, but shiftes to larger values as v tends 
to 218 GHz. Note, however, that at this frequency the tSZ signal 
drops to zero. 



Correlation Coefficient R 
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Correlation coefficients 7ti = 
)l\JCf MB C\ sz for clusters in the ACDM (thick 



line) and SCDM (thin line) scenarios. 



ACDM and SCDM cosmological models (thick and thin 
solid lines, respectively), (figure (0J). The ISW is the 
cause of the coupling of CMB anisotropies with tSZ 
signal in the ACDM case. This causes a cross-correlation 
with the total CMB signal of about a 20% at I ~ 10, 
which drops at higher multipoles since the ISW signal 
decreases rapidly with increasing I. For the SCDM model, 
we obtain l min ~ 5 for rjtsz{z ~ 0.5) ~ 6730 Mpc and 
rjo ~ 8300 Mpc; and l ma x ~ 8, which would explain the 
low level of correlation in this case (less than a few percent). 



3.2 Reionization and Resonant Scattering of 

CMB Photons on Ions, Atoms and Molecules 
of Heavy Elements 

Both scattering on free electrons during reionization and 
resonant scattering associated to any type of transition in 
heavy species contribute with some optical depth for the 
CMB photons. In the first case, the optical depth is gen- 
erated by the Thompson scattering occuring between CMB 
photons and free electrons, and, hence, is frequency indepen- 
dent. This situation changes for resonant transitions, pro- 
vided that CMB photons scatter the line only if their fre- 
quency is close enough to the resonant frequency. Apart from 
this distinction, the effect of both phenomena on the CMB 
power spectrum is identical, so we shall restrict our analysis 
on the case of resonant scattering, (which, by its spectral 
peculiarity, can be separated from the intrinsic CMB tem- 
perature fluctuations). Hence, we refer to Basu, Hernandez- 
Monteagudo & Sunyaev (2004), (hereafter BHMS) where 
this effect is utilized to discuss constraints on the abun- 
dances of heavy species at redshifts 0.1 < z < 30. 

If we denote by r rs the homogeneous (i.e. position inde- 
pendent 7 ) optical depth associated to resonant scattering, 
we can write that the change induced by the resonant tran- 
sition on the temperature field is given by: 



T r . 



ib e 



■To. 



(27) 



where T rs is the temperature angular fluctuation field at the 
time of resonant scattering, T cm (, is the intrinsic CMB field 
generated at the LSS and T gen are the new temperature 
fluctuations generated by the resonantly scattering species. 
If we now take the limit r rs <C 1, the last equation becomes 



T rs — (1 7Y S ) T C mb ~b T r s Tg en -f" 0[T rs 



(28) 



with Tg™n the coefficient of the linear term in the expansion 
of Tg e n in terms of r rs . In Fourier space, this translates into: 

\ , lin 

&(Z,m), rs — V-L 7~rs) a (l,m), crab T Trs <2(/ im ) gen 

+ 0[r 2 rs ], (29) 

with a(; !m )'s denoting Fourier multipoles. If we now define 



l (l.m) 

SCi = (|o(j )T1 
find that 



,fV<|o(J, 



0, 



%b \ ) , it is straightforward to 



SCi = r rs 2 {{Tie [o(j, m ), 



cmfe X (ttf^ m ^ gen ) J ) 
(l^(i,m), cmb 



2 >) + 0[T 2 rs 



(30) 



As shown in detail in the Appendix A of BHMS, the first 
term accounts for the correlation between fluctuations gener- 
ated during recombination and those generated in the epoch 
of resonant scattering, whereas the second (autocorrelation) 
term expresses the blurring of the intrinsic anisotropies in- 
duced in the LSS due to the resonant scattering at lower 
redshift; from now this term will be referred to as the blur- 
ring term. Note that it is merely proportional to the intrin- 
sic CMB power spectrum at the resonant scattering epoch, 
and hence, as long as resonant scattering takes place after 



7 In the optically thin limit, r rs <g 1, one can relax the approx- 
imation on homogeneity by assuming that the scales at which 
T rs varies are smaller than the scales under study, for which an 
average integrated optical depth is effectively working. 
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recombination, the shape of this blurring term will be iden- 
tical to the primordial CMB power spectrum generated at 
decoupling, and thus redshift independent. For the reasons 
outlined at the end of Section 2, the correlation term is only 
of relevance at the very low I range of multipoles, in which 
newly generated anisotropies overcome the blurring of origi- 
nal temperature fluctuations and introduces new anisotropy 
power, (see again Appendix A of BHMS). This occurs for 
both ACDM and SCDM cosmological models, since the In- 
tegrated Sachs- Wolfe effect (hereafter ISW) has no effect 
here provided that, in adiabatic A models, it becomes im- 
portant only at very low redshift, during the A term domi- 
nance, whereas for an Eins tein-de Sitter Universe i t vanishes 
in the linear regime, (e.g.. lHu fc Sueivamal dl995T) 'l. Recall- 
ing that r/rec ~ 300 Mpc, r/ r3 (z — 25) ~ 2811 Mpc, and that 
770 — 14000 Mpc, one finds that l m in ~ 1 (no drop of the 
cross correlation expected at low multipoles) and l ma x ~ 30, 
at which we would expect having some decrease in the am- 
plitude and/or change of sign in the cross correlation. Figure 
shows the actual computation of the terms in eq. 13011 : all 
curves have been computed for r rs — 10~ 3 and rescaled to 
r rs = 1, so the actual measurement that our method would 
provide is then given by the diamonds line times r rs , (which, 
for small enough r rs , is below the cosmic variance limit). As 
in BHMS, the resonant lines have been modelled by a gaus- 
sian centered on the conformal time (r) rs ) corresponding to 
the redshift considered in each case, and with a a equal to 
one percent of rj rs . Solid lines gives the blurring of the orig- 
inal power spectrum, and the dashed line accounts for the 
cross-correlation. Note that we are plotting absolute values, 
and that only at low multipoles the first term is positive 
and greater in amplitude than the blurring term. For higher 
multipoles, the correlation term can be neglected and one is 
left with the simple autocorrelation term: 

SC t ~ -2r rs Cr"- (31) 

This ^-dependence for the SCi's is generic for any source of 
localised optical depth for the CMB photons. This drop in 
the SCi's (and change of sign at some l zero , see low / range 
in figure JSJ) are a direct consequence of the correlation of 
fluctuations at r) rec ,rj rs , and provide a test for the origin 
of the <5CVs, just as in the case of the tSZ effect addressed 
above. 

We remark that these SCi's are measurable only if the 
CMB is being observed at two different frequencies; one cor- 
responding to the resonant scattering at r/ r3 , and another 
one in which such resonant scattering can be neglected. Note 
that there is no place for this situation in the case electron 
scattering during reionization, since Thompson scattering on 
free electrons is frequency independent. We are implicitely 
assuming that the instrument is sensitive to the amplitudes 
of the SCi's: in BHMS we showed that the current detec- 
tor technology (present in experiments like WMAP, ACT 
or Planck) should already allow to set strong limits in the 
abundance of resonant species during the epoch of reioniza- 
tion. 

3.3 Emission in Fine Structure Lines of C, N, O 
in Haloes 

BHMS studied the effect of resonant scattering of 
CMB photons in fine structure transitions associated 
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Figure 5. Angular power spectrum arising as a consequence 
of resonant scattering on a line placed at the end of the Dark 
Ages, (2 = 25). We are plotting the two terms contributing to 
the total SCfs (diamonds): the absoption term is displayed in 
thick solid line, and is merely proportional to the intrinsic CMB 
power spectrum. It is negative, and hence the total SCi's cross 
zero only when the Doppler-induced generation (dashed line) be- 
comes relevant at low multipoles. Only absolute values are shown. 
The correct amplitude of this effect is obtained after multiplying 
these curves by T rs - 



to metals and ions. They found that very overdense 
regions (S 10 4 ) should emit in these lines via 

collis i onal excitations, iSuginohara. Suginohar a, fc SpereeJ 
ll999UVarshalovich. Khersonskii. fc Sunvaevlll98lll . The ex- 
pected amplitude of this signal is relatively small, while its 
spectral dependence is very different from that of the CMB. 
On the other hand, it also depends on the star formation 
history in haloes whose large scale distribution should trace 
the general density fluctuation field. For these reasons, one 
can consider the application of the correlation method in 
this case as well. The main difference to the scenario stud- 
ied by BHMS is that, in this occasion, the scattering in the 
lines is almost negligible, and hence, no blurring of original 
CMB anisotropies should be expected. Hence, there will be 
no further suppresion of the CMB power spectrum at high 
multipoles, but only extra power in the large angular scale 
range. This is motivation of an upcoming paper where both 
the linear and quadratic terms are taken into account. 

3.4 Extragalactic Foregrounds 

In this subsection we address possible effects that well- 
known physical processes (such as free-free emission, dust 
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emission in the IGM or inside galaxies and synchrotron 
emission in extragalactic radio sources) have on our method. 
In the case of extragalactic foregrounds, it is clear that if 
they are produced in haloes, they should trace the overall 
mass distribution in the very large scales, just as in our 
study of tSZ signal induced by clusters of galaxies. For this 
reason, one could think of applying this method on them, 
expecting to find a similar shape for the correlation term 
at large angular scales as the one found for tSZ clusters. 
This raises the question whether these foregrounds could 
mutually contaminate or bias the correlation estimates. 
Since the method proposed here is based on the frequency 
dependence of the signal under study, proper frequency 
coverage should allow to identify and separate each compo- 
nent as long as spectral signatures are distinct enough. 

It is obvious that if the sources of these signals are lo- 
cated in our Galaxy, one would not expect any type of cor- 
relation between them and the original density perturbation 
field, leading to no linear (foe e) term. 



4 DISCUSSION AND CONCLUSIONS 

The amplitude of the cross correlation depends essentially 
on the conformal distance separating the signal sources, 
rather than the particular k projection of sources of different 
origin. The closer the sources of the signals are, the higher 
the correlation becomes. At this respect, the presence of a A 
term generating an ISW signal is of crucial importance for 
those effects g enerated in ou r neighborhood, (particularly 
the tSZ effect. ICooravl feOOlb '). Consistently with the ISW 
contribution to the total CMB signal (around 20 /iK with 
respect the total ~ 110 /iK of the CMB), the correlations 
in a ACDM universe show typical values of 10-20 %, with 
remarkably lower values in the SCDM scenario. In an 
Einstein-de Sitter universe, the correlation drops to a few 
percent, and the enhancement of the weak signal is rather 
far from being relevant. The situation changes remarkably 
in the case of resonant scattering at high redshift. In this 
situation, the correlation coefficient is practically unity 
for the low multipoles, since, as shown in BHMS, arises 
as a consequence of the monopole and Doppler terms of 
the CMB, and the contribution of the ISW component is 
negligible. 

Although the galactic contamination is thought to 
be more important in the large angular scales where 
these correlations show up, it is also expected that space 
experiments achieve their best sensitivities in the big 
angular scales. In the case that the signal is of extragalactic 
origin, the cross term will always show up together with 
the squared term, although both terms have, in general, 
different frequency 112MB and /-dependence. This should also 
help in distinguishing between them, specially in the case 
of the tSZ effect, for which a peculiar pattern of the <5CVs 
versus / and v has been predicted. In the high multipole 
range, frequency dependent scattering such as resonant 
scattering introduce a measurable blurring of original CMB 
temperature fluctuations generated during recombination. 
Since it merely consists in an autocorrelation of CMB 
anisotropies, this blurring term has the same 1-dependence 



as the original CMB power spectrum. 

The method proposed here can also be applied in the 
study of the cross correlation of CMB temperature fluctua- 
tions with the radio background. In the low frequency range, 
new instruments like the Low Frequency Array (LOFAR) or 
the Square Kilometer Array (SKA) will measure the radio 
background. This is mainly due to radio galaxies present 
in the redshift range z £ [0,4], and its fluctuations are 
expected to be of much higher amplitude than those of the 
CMB. However, due to the fact that the radio background 
is generated by radio galaxies tracing the universal density 
fluctuation field, one can think of applying this method in 
order to enhance the CMB component at these frequencies. 
When doing this, one must keep in mind that there is 
emission at 21 cm coming from neutral hydrogen during 
the D ark Ages, (z ~ (30, 100). iMadau. Meiksin fc Reel 
lll997t) ) which should fall in this frequency range and which 
is showing also some degree of correlation with the CMB. 
However, according to the arguments given in Section 2, 
most of the correlation will be due to the coupling of the 
ISW effect with the radio galaxy distribution at low and 
moderate redshifts. 

Similar arguments can be applied when studying the 
857 GHz band of Planck's HFI, since we can expect that 
this method should be able to unveil the distribution of 
extragalactic dust and its imprint on the CMB. In other 
words, by means of the ISW the CMB has become a tool 
which permits performing independent tests at different 
frequencies on the large scale distribution of matter. The 
main two caveats to have present are the possibility of 
having some signal generated during reionization, at very 
high redshift, which could be introducing some extra 
correlation, and the presence of galactic foregrounds, whose 
residuals might invalidate these analyses in the very low 
multipoles. 

In this paper, we have addressed the issue of correlated 
signals in the context of CMB. We have shown that, 
in the case in which two signals have different spectral 
dependence, the presence of correlations between both 
can be used in order to enhance the weak signal with 
respect the dominant one. Assuming that the correlation 
between signals is caused by the the cosmological den- 
sity perturbation field, we have found at which angular 
range such correlation might be relevant. This depends 
essentially on two different scales: the distance separating 
the events generating the signals under consideration, 
and their distance to the observer. In a ACDM universe, 
these cross terms dominate at the large angular scales, 
and hence characterize our predictions of the power spec- 
tra associated to the weak signals in the low multipole range. 
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